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Abstract. Galaxy clusters are valuable cosmological probes. However, cluster mass estimates 
rely on observable quantities that are affected by complicated baryonic physics in the intracluster 
medium (ICM), including feedback from active galactic nuclei (AGN). Cosmological simulations 
have started to include AGN feedback using subgrid models. In order to make robust predictions, 
the systematics of different implementations and parametrizations need to be understood. We have 
developed an AGN subgrid model in FLASH that supports a few different black hole accretion 
models and feedback models. We use this model to study the effect of AGN on X-ray cluster 
observables and its dependence on model variations. 
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INTRODUCTION 

Clusters of galaxies are useful probes of cosmological parameters, provided that their 
masses can be determined accurately from observables such as X-ray luminosity or 
X-ray temperature. However, it is still a challenge for current theoretical models to 
reproduce all the observed ICM properties. Feedback from active galactic nuclei (AGN) 
has been commonly invoked as the necessary missing piece. Since AGN feedback 
involves a wide range of scales, from the accretion disk on AU scales to clusters on 
Mpc scales, only recently has it been included in cosmological simulations using subgrid 
modeling. 

Because there are a variety of implementations and parametrizations of AGN subgrid 
models in the literature, the aim of this study is to understand the current theoretical 
uncertainties in predicting the ICM properties. More specifically, we compare two black 
hole (BH) accretion models and two feedback models that are commonly used: bub- 
ble feedback [1] and jet feedback |@]. We will describe the analytical and numerical 
approaches in the following section. In Section 3, we will discuss the effect of AGN 
feedback on cluster observables with these model variations. 



METHODOLOGY 



We performed hydrodynamic simulations with radiative cooling and AGN feedback 
with an isolated cluster sitting in a 2048 kpc box using FLASH 2.5 yfl. The region 
surrounding a central 10 7 ft M© black hole is refined to the maximum resolution of 4 
kpc. Following [4], the cluster gas is initialized assuming a modified NFW [5] profile and 
is in hydrostatic equilibrium in an NFW gravitational potential. The cluster has a mass 
of 10 14 ft M©, concentration of 6.5, cool-core radius of 24.5 kpc, and gas fraction of 
0.12. Radiative cooling is computed using ((J assuming zero metallicity. The Hubble 
constant h = 0.7 is used. 

We relate the BH accretion rate to the Bondi-Hoyle-Lyttleton |0] accretion rate: 

M B h « M Bo ndi = 4nG 2 Ml n p/c 3 sl (1) 

where A/bh is the BH mass, and p and c s are the gas density and sound speed, respec- 
tively. The proportionality reflects the fact that cosmological simulations usually do not 
have sufficient resolution to resolve the Bondi radius, and hence using only the Bondi 
accretion rate would underestimate the actual BH accretion rate. There are a number of 
ways to link these two quantities. In this study we will present results using the constant- 
ce model, A/bh = ocMBondi AID' an ^ the density-dependent /3 model [8|]: 

Mbh = (n H /n* H ) p M Bondu (2) 

where n* H = 0. 1 cm~ 3 . We impose an upper limit on the accretion rate corresponding to 
the Eddington rate, M Edd = (4^GMbhw p ) /(gfOrc), where m p is the mass of the proton, 
Oj is the Thompson cross-section and £f is the radiative efficiency. 

The feedback from the AGN to the surrounding gas is linked to the BH accretion rate. 
For bubble feedback, as in [1], the bubbles are only formed when the BH mass increases 
by a fraction 5bh since the last bubble formation. When a bubble is formed, only thermal 
energy is injected: 

E = £ m £ f c 2 5M BH , (3) 

where £ m is the efficiency of mechanical heating, and 5A/bh = ^bhA^bh is the increase 
in BH mass since the last bubble was formed. The injected energy is evenly distributed 
within a sphere of radius, 

where the scaling parameter values Rq = 30 h ! kpc, Eq = 5 x 10 erg, and po = 
10 6 /? 2 M Q kpc 3 are motivated by observed bubble sizes. The bubbles are centered on 
the black hole. The default bubble parameters are a = 350, /3 = 2, £f = 0.1, £ m = 1, and 
5bh = 0.001. 

For jet feedback, we adopt the model in |2|], where the injection rates of the mass, 
momentum, and energy are given by 

M = TfM BH m (5) 
P = v^Mbhc^, (6) 
E = £ f M BH c 2 m (7) 



where r\ is the mass loading factor, and the function \P determines the spatial extent of 
the jet: 

Y(x) = — ^exp \-^-\ JL (8) 

The jet is aligned with the z-axis and the feedback is applied to regions with |z| < h eq 
and r < 2.6r e j. The default values of the jet parameters are r\ = 100, £f = 0.1, r e j =3.2 
kpc, and /z e j = 2.5 kpc. 



RESULTS 



Figure 1 shows the growth of the black hole and its accretion history for different 
accretion and feedback models. We can see that the overall behavior is similar: the black 
hole accretion rate (BHAR) is regulated after t ~ 0.5 Gyr, which leads to a converged 
BH final mass that is in general consistent with observed BH masses in 
clusters. However, the paths of BH growth are different, namely, jet feedback results 
in a more gradual increase in BH mass, while for the bubble feedback the black hole 
grows more rapidly initially and reaches a quasi-static state earlier than in the jet model. 
This difference is mainly driven by the behavior of the accretion rate at the early stage 
(t < 1 Gyr). Since the bubbles are only allowed to form when the BH mass exceeds the 
threshold fraction 5bu, the black hole is able to accrete more before feedback begins to 
dominate. In contrast, the continuous jet feedback slows down the accretion from the 
beginning and thus yields a more mild BH accretion and mass growth. Similarly, the 
black hole in the /3 model grows a little more slowly than in the a model due to the 
difference in the initial estimate of the accretion rate. 
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FIGURE 1. Evolution of BH mass (left) and BH accretion rate in units of the Eddington rate (right) for 
different accretion and feedback models. 

In Fig. 2 we show the evolution of cluster X-ray observables. We can see that different 
accretion and feedback models actually can have quite different results. We find that in 
terms of heating the ICM, the a model with bubble feedback is most effective and the 
jet feedback is least effective. One possible reason is that since the amount of injected 
energy is proportional to the BH mass, less feedback energy is available for the jet model, 
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FIGURE 2. Evolution of X-ray luminosity (left) and spectroscopic-like temperature (right) for different 
accretion and feedback models. 



in which the black hole grows more slowly. Moreover, only a portion of the feedback 
energy is turned into thermal energy in the jet model, while in the current setup for 
bubble feedback, all the energy is used to heat the cluster. Using a more realistic e m may 
reduce the discrepancy. 

In conclusion, using a set of isolated cluster simulations we find a non-negligible 
difference in predicting the growth of black hole and influence on cluster observables 
using different subgrid AGN models. Model uncertainties should be carefully examined 
in future cosmological simulations with AGN feedback. 

ACKNOWLEDGMENTS 

The authors acknowledge support under a Presidential Early Career Award from the U.S. 
Department of Energy, Lawrence Livermore National Laboratory (contract B532720). 
Additional support was provided by NASA Headquarters under the NASA Earth and 
Space Science Fellowship Program (NNX08AZ02H), a DOE Computational Science 
Graduate Fellowship (DE-FG02-97ER25308), and the National Center for Supercom- 
puting Applications. The software used in this work was in part developed by the DOE- 
supported ASC / Alliance Center for Astrophysical Thermonuclear Flashes at the Uni- 
versity of Chicago. 

REFERENCES 

1. D. Sijacki, V. Springel, T. Di Matteo, and L. Hernquist, MNRAS 380, 877-900 (2007). 

2. A. Cattaneo, and R. Teyssier, MNRAS 376, 1547-1556 (2007). 

3. B. Fryxell, et al., Astrophys. J. Supp. 131, 273-334 (2000). 

4. D. Sijacki, and V. Springel, MNRAS 366, 397-416 (2006). 

5. J. F. Navarro, C. S. Frenk, and S. D. M. White, MNRAS 275, 720-740 (1995). 

6. R. S. Sutherland, and M. A. Dopita, Astrophy. J. Supp. 88, 253-327 (1993). 

7. H. Bondi, MNRAS 276, 663 (1952). 

8. C. M. Booth, and J. Schaye, MNRAS 398, 53-74 (2009). 



10.00 




0.01 



